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The interaction betweenT-cell receptors (TCRs) and peptide epitopes is highly degenerate: 
aTCR is capable of interacting productively with a wide range of different peptide ligands, 
involving not only cross-reactivity proper (similar epitopes elicit strong responses), but also 
polyspecificity (ligands with distinct physicochemical properties are capable of interacting 
with theTCR). Degeneracy does not gainsay the fact thatTCR recognition is fundamen- 
tally specific: for the vast majority of ligands, the functional sensitivity of a given TCR is 
virtually null whereas this TCR has an appreciable functional sensitivity only for a minute 
fraction of all possible ligands. Degeneracy can be described mathematically as the proba- 
bility that the functional sensitivity, of a given TCR to a randomly selected ligand, exceeds 
a set value. Variation of this value generates a statistical distribution that characterizes TCR 
degeneracy. This distribution can be modeled on the basis of a Gaussian distribution for 
the TCR/ligand dissociation energy. The kinetics of theTCR and the MHCI molecule can 
be used to transform this underlying Gaussian distribution into the observed distribution 
of functional sensitivity values. In the present paper, the model is extended by accounting 
explicitly for the kinetics of the interaction between the co-receptor and the MHCI mole- 
cule. We show thatT-cells can modulate the level of degeneracy by varying the density of 
co-receptors on the cell surface. This could allow for an analog of avidity maturation during 
incipient T-cell responses. 
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INTRODUCTION 

Thymus-derived lymphocytes (T-cells) recognize peptide antigens 
via antigen-specific receptors (TCRs); in particular, CD8 + cyto- 
toxic T lymphocytes (CTLs) recognize short peptides presented 
by major histocompatibility complex (MHC) class I molecules 
(1). Estimates of the human TCR diversity suggest that there are 
~10 8 different antigen receptors in the naive T-cell pool (2), which 
raises the question how such a limited TCR repertoire can provide 
effective immunity to perhaps over 10 15 distinct pMHCs (2). The 
discrepancy suggests that even while TCRs are highly specific, a 
considerable degree of degeneracy remains (we use degeneracy as 
a term of art to cover both polyspecificity and cross-reactivity). 
The central importance of degeneracy was first pointed out by 
Mason (3) and later confirmed by others such as (4). Experi- 
mental and mathematical studies confirm that TCR recognition 
is highly degenerate: a single TCR may be able to recognize, at 
physiologically relevant degrees of functional sensitivity, over one 
million different peptides in the context of a single MHCI molecule 
(2, 5), an estimate that takes into account the binding specificity 
of the MHC molecule, but not the additional selection imposed 
by the stringency requirements of peptide cleavage and process- 
ing in the presentation pathway. The latter constitutes an epitope 
diversity filter that is instrumental in regulating the immuno- 
visibility of salient epitopes (6). The issue of whether there exists 
an optimal level of TCR repertoire diversity was reviewed by 



Nikolich-Zuglich et al. (7), and various authors have reviewed 
the functional repercussions at the systems level (8-10). 

The interaction between TCR and pMHCI ligand can be mod- 
ulated by the co-receptor CD8 in several ways: (i) promoting the 
association of TCR and pMHCI; (ii) stabilizing the TCR/pMHCI 
interaction; and (iii) enhancing the rate at which the TCR/CD3 
complex attains signaling status by association of TCR/CD3 with 
protein tyrosine kinases such as p56 lck and adaptor molecules such 
as LAT and LIME (5, 1 1-16). The first of these three mechanisms 
modifies the affinity of the TCR/pMHCI interaction, whereas the 
third alters the time it takes for an engaged TCR/CD3 com- 
plex to attain full signalosome status. In particular, CD8 can 
enhance the TCR/pMHCI association rate by 50%, and reduce 
the TCR/pMHCI dissociation rate by at least 50% (16, 17), and 
CD8 modulates the rate of immune receptor tyrosine-based acti- 
vation motif (ITAM) phosphorylation, by recruiting TCR/pMHCI 
complexes to membrane micro-domains at a rate which depends 
on the affinity of pMHCI/CD8 binding (16). 

These findings suggest that CD8 not only controls degener- 
acy, but also differentially regulates functional sensitivity, that is, 
the T-cell can increase its sensitivity for one ligand, while reduc- 
ing it for others. By varying the level of CD8 expression, the 
T-cell can increase its sensitivity to the disease-associated antigen, 
while at the same time decreasing its sensitivity to antigens asso- 
ciated with healthy conditions. This novel mode of co-receptor 
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FIGURE 1 |TCR/pMHCI/CD8 interactions represented in the model 

(A) Kinetics of TCR/pMHCI/CD8 complex formation. TCRs and CD8s are 
located on theT-cell side and pMHCIs on the APC side. Four possible states 
are shown: in (I) and (II), pMHCI is not bound to CD8 whereas in (III) and 
(IV), pMHCI is bound to CD8; in (I) and (III), pMHCI is not bound toTCR/CD3 
whereas in (II) and (IV), pMHCI is bound to TCR/CD3. Transition rates 
between these four states are indicated with forward rates labeled as X, 
and backward rates labeled as X_, (where / = 1, 2, 3, 4). The white circles in 
the tails of theTCR/CD3 complex represent the ITAMs, which are present 
on the cytoplasmic tails of the y, 8, the two e, and the two i; polypeptide 
chains of the CD3 complex; only one such tail is indicated here, for clarity. 
There are n such ITAMs, with only three shown here for the sake of clarity. 
This part of the kinetic scheme is assumed to be in equilibrium. (B)The 
TCR/pMHCI bound states (II) and (IV) are further subdivided into n states 
each depending on the number of phosphorylated ITAMs, resulting in a 
Markov chain. The triggered states are the right-most ones where all n 
ITAMs have been phosphorylated (indicated as filled-in circles). At each 
state a reversal (at rate X_, or X_J to the unbound states (I) or (III), 
respectively, can occur. Transitions between the sequences occur at rates 
X 2 and X_2- Left-to-right transitions within the sequences correspond to 
ITAM phosphorylation by kinases such as p56 lck and ZAP-70, at rate 
nX = n/T R in the sequence with CD8 unbound and at rate nX* = n/T* in the 
chain with CD8 bound. The Markov chain serves to calculate the TCR 
triggering probability and is not assumed to be in equilibrium. 



action could be critical in ensuring that the TCR repertoire retains 
the ability to respond to antigenic challenges, while avoiding 
autoimmunity. 

T-cell antigen recognition can be expressed in terms of its func- 
tional sensitivity (18). One of the main determinants of functional 
sensitivity is the rate at which a single agonist copy elicits TCR 
triggering. Functional sensitivity depends on bio-molecular para- 
meters such as the TCR/pMHCI on-rate and off-rate. The mole- 
cular kinetics at the T-cell:antigen-presenting cell (APC) interface 
determine this relationship. This kinetic theory resolves the long- 
standing controversy over whether T-cell activation is governed 
by affinity or off-rate (cf. (19, 20)); the theory shows that both 
parameters play a role, but in the so-called MHC-limited regime, 
the off-rate is the main governing factor. However, the on-rate 
and the off-rate together determine whether or not the kinetics is 
MHC-limited (21-23). 

The aim of the present study is to explore the kinetic basis 
of the role of CD8 in regulating degeneracy and functional sen- 
sitivity. Our model generalizes the classical kinetic proofreading 
model, as proposed by McKeithan (24), and later modified by 
others (21, 22, 25-27). The kinetic proofreading model assumes 
that the TCR needs to remain bound to a pMHC in order for the 
TCR/CD3 complex to become fully activated, via several modifi- 
cation steps such as phosphorylation of several tyrosine residues 
on the TCR complex, and recruitment and subsequent activation 
of ZAP-70 (28). 

The statistics of TCR degeneracy is modeled by relating the 
TCR/pMHCI mean interaction time with the dissociation energy 
according to Arrhenius theory. We show that by varying the total 
co-receptor density and key kinetic parameters, the T-cell can 
modulate the level of degeneracy. Furthermore, we compare our 
results with experimental data for HLA A*0201 mutants with 
altered binding affinity for CD8. 

THEORY: TCR/pMHCI/CD8 KINETICS AND TCR DEGENERACY 

We assume (i) that the TCR/CD3 complex on the T-cell sur- 
face becomes triggered (achieves signalosome status) during an 
interaction with a pMHCI ligand if it undergoes n ITAM phos- 
phorylations, where n is a positive integer; and (ii) that the 
presence of a CD8 molecule bound to the pMHCI ligand affects 
the TCR/pMHCI association and dissociation rates, as well as 
the rate at which the TCR/CD3 complex progresses through 
the phosphorylation sequence. The co-receptor CD8 is associ- 
ated with TCR/CD3 with tyrosine kinases such as p56 lck , which 
phosphorylates tyrosine residues within the ITAMs. 

To determine the kinetics of the TCR/pMHCI/CD8 inter- 
actions, consider the four binding states of a pMHCI mole- 
cule, labeled (I-IV) in Figure 1A. The states in which the TCR 
and the pMHCI ligand are bound can be further subdivided 
into different states according to the number of phosphorylated 
ITAMs, as shown in Figure IB. The transition rate between 
these states is h/Tr, where Tr is the average time required to go 
through the entire sequence. The notation Tr emphasizes that 
this quantity is the average time for the TCR/CD3 complex to 
become triggered. We call this the TCR triggering threshold. Typ- 
ical values of TCR triggering threshold Tr are in the range of 
5-15 s (29). 



KINETIC EQUATIONS 

The model describes the kinetics of the interactions between the 
TCR, pMHCI molecules, and co-receptor CD8 in the contact area 
between a T-cell and an APC. This area is occupied by TCRs and 
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CD8s on the T-cell side and by pMHCIs on the APC side, whereby 
CD8 binds pMHCI at a distinct site of the TCR/pMHCI complex. 
ITAMs are located on the ^-chains of the CD3 complex associated 
with the TCR and are represented schematically in Figure IB. 

Interactions between TCR, pMHCI, and CD8 

Key quantities in the model (summarized in Table 1) are the den- 
sities of free receptors on the T-cell and free ligands on the APC, 
which can also form complexes with the receptors. The surface 
densities of the TCR, CD8, and pMHCI molecules are subject to 
the following conservation laws: 



Table 1 | Model parameters and variables. 



M 


Free pMHCI density 


Mr 


TCR/pMHCI density without CD8 bound 


M x 


pMHCI/CD8 density without TCR bound 


M X R 


TCR/pMHCI/CD8 density 


X 


Free CD8 density 


R 


Free TCR density 


M T 


Total pMHCI density 


X T 


Total CD8 density 


Rt 


Total TCR density 



Mj M + Mr + Mx + Mxr , ( 1 ) These results can also be expressed in terms of standard affinity 

X T = X + M x + M XR , (2) constants: 

R T = R + Mr + M X r , (3) 



M T R 

M R = 



which state that for each species, the total density must equal the K x + R + XK X /K 3 + RX/K 2 

sum of all bound forms plus the free density. We define pseudo- ^ _ MjX 

unimolecular association rates X/(i = 1, . . ., 4) as follows (see K 3 + X + RK 3 /K\ + RX/K4 

Figure 1A): M T XR 

M X r ~ 



(13) 



K X K 2 + RK 2 +XK4 + RX 
X, = A\R, X 2 = A 2 X, (4) 

X 3 = A 3 X X4 = A4.R (5) ^he co-receptor CD8 modulates the rate of TCR triggering. 

Three major modulatory functions of the co-receptor have been 
where A,-A 4 are the two-dimensional association rates documented: modulation of TCR/pMHCI on-rate, TCR/pMHCI 
(cmV 1 ) for TCR/pMHCI or CD8/pMHCI binding. Two- off-rate, and of the ITAM phosphorylation rate. These effects can 
dimensional dissociation constants (cm" 2 ) are defined as follows: be represented by dimensionless multipliers. 

(i) enhanced TCR/pMHCI on-rate: 
A 4 = YonAi where y on > 1; 

(ii) reduced TCR/pMHCI off-rate: 
^•-4 = YofA-i where 0 < y off < 1; 



K! = -±R, 

M 


(6) 


K 2 = —±X, 


(7) 


^3 = -X, 

^3 


(8) 


K4 


(9) 



In order for a system of reactions to be in thermal equilibrium, 



(iii) increased phosphorylation rate, which is equivalent to a 
reduced TCR triggering threshold T| 



each individual reaction must be at equilibrium (the principle of x = ^ where yR < : and X = 1 and = 1 . (14) 

detailed balance): R J « 

XiM — X iMr ~K 2 Mr — X 2 Mxr It is sometimes convenient to combine the on-rate and off-rate 

(10) effects into a single coefficient, as follows: 

X 3 M = X_ 3 M X , X 4 M X = X-iMxR, 

Yoff K4 K 2 

from which it follows that K 1 K 2 = K 3 K i . Combining this with Ykin = — = — = —, where 0 < y kin < 1. (15) 

Yon -*^3 

the conservation laws with we obtain: 

Mj Mj\\ From equation (15), we have < y Q g. We rewrite Mr, Mx, 

M = , Mr = , and Mxr, as given by equation (13), in terms of Yki n : 

C C ^ (11) 

M T X 3 M T \i\ 2 r/k 

x = OvZ7' XR = cxTTxT^' M * = M r 



1 + R/Ki + X/Ki + RX/ (K1K3 Yki„) 



X/K } 

where Mx = Mj , (16) 

1 + R/Ki + X/K 3 + RX/ (^1^3Ykin) 

Xi X 3 XiX 2 , , RX/(KiK 3 ) 

C= 1 + — - + -^+ - \ l . (12) M XR = M T - 



X-i X_ 3 X_,X_ 2 > y kin( i + X /K 3 + R/K l )+ XK 3 /(RKi) 
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To non-dimensionalize this system, we introduce the following 
dimensionless quantities: 



k = 



X 



R 

K~i 



XT = 



Xj 



r T = 



R T 



mj = 



Mt 



(17) 



It follows from equation (17) that MjlKj, = Krm. By equa- 
tions (16) and (17), we have the following non-linear system of 
equations that determines the surface densities of free and bound 
receptors of all species: 



XT = X+KtTlT- 



xr 



1 + x + r + xr/y kin y kin (l + x + r) + xr 

r xr 

r T = r + fflj- h raj . 

1 + x + r + xr/y kin y kin (l + x + r) + xr 

(18) 

This system is readily solved numerically for x and r, given the 
parameters y k { n , k, rj, and mj. 



The TCR triggering rate 

The functional sensitivity of the TCR is represented in the present 
model as the rate at which TCR/CD3 complexes attain signalo- 
some status. To calculate this TCR triggering rate, which we shall 
denote by W, consider the Markov chain as depicted in Figure IB, 
which is best thought of as a system of two coupled linear Markov 
chains. The TCR/CD3 complex has to undergo n phosphoryla- 
tions for the TCR to be triggered. When the co-receptor is not 
bound to the TCR/pMHCI complex, individual phosphorylation 
steps proceed at rate «/Tr = rik whereas when the co-receptor is 
engaged, this rate is n/T^ = rik* > rik. The factor n arises simply 
as a matter of scaling, so that the average time to progress through 
the chain of phosphorylations equals Tr. 

In reality, signalosome formation involves several other types 
of event besides ITAM phosphorylation, such as binding of ZAP- 
70, engagement of LAT, and so on. To avoid cumbersome notation 
we shall formulate the model as if ITAM phosphorylation were the 
only type of event; the essential theory is not materially affected by 
this simplification. We do not assume an equilibrium state for the 
Markov chain: the complex starts at zero phosphorylations at the 
beginning of every encounter with a pMHCI ligand and proceeds 
forward stochastically. 

The pMHCI/TCR/CD3 complex may not attain the nth state 
(triggered state), but instead the pMHCI/TCR engagement may 
terminate, which happens at rate X_i when the co-receptor is not 
bound and at rate X_4 < X_i when the co-receptor is engaged. 
We assume that upon TCR/pMHCI dissociation the CD3 com- 
plex reverts to the basal state of zero ITAM phosphorylations 
sufficiently rapidly that the TCR/CD3 complex will be in this com- 
pletely unphosphorylated state when the next encounter with a 
pMHCI molecule occurs. Essentially, this means that the CD3 
complex is more susceptible to the action of phosphorylases 
and/or less susceptible to the action of kinases when the TCR 
is not engaged. A mechanistic explanation underpinning this 
assumption lies outside the scope of the present model. 



The probability that the TCR/CD3 complex will undergo 
another ITAM phosphorylation is given by 



rik 



X-i + n\ + X 2 
in the case where the co-receptor CD8 is not engaged and by 



(19) 



X_ 4 + n\* + X_ 2 



(20) 



when CD8 is engaged. The following expression gives the proba- 
bility that the system switches from the chain with CD8 unbound 
to the chain with CD8 bound: 



1 = 



X. 



X 2 

-i + nk + X 2 



(21) 



whereas the system will switch from CD8-unbound to CD8-bound 
with probability 



X-2 



X_ 4 + n\* + X_ 



(22) 



TCR triggering requires completion of all steps before the 
TCR/pMHCI complex comes apart. We shall find an expression 
for the probability that the CD3 complex will ultimately attain 
completion when starting from i completed steps. Let P° denote 
this probability for the case with CD8 unbound and P* for the case 
with CD8 bound. The TCR triggering probability is then found 
as P° if CD8 is unbound when pMHCI docks the TCR, and P^ if 
CD8 is bound. These triggering probabilities allow us to calculate 
the TCR triggering rate W. In particular, W can be expressed as the 
rate at which TCR/pMHCI complexes dissociate, times the prob- 
ability that whenever a given TCR/pMHCI docking commences, 
the CD3 complex is ultimately triggered (the probabilities Pq and 
Pq )• In formula, this statement is represented as follows: 



(23) 



which, by the principle of detailed balance, can be rewritten in 
terms of dissociation rates, as follows: 



W = R\- 1 W a 0 + M XR \- i P* 0 . 



(24) 



The law of total probability yields the following system of 
coupled difference equations: 



i-l 
i-l J 



1 



1^ 



b 11 a 



rpn 

Ln J 



which can be solved to give 



(1 





pop*" 


n 


' i " 








p*p0 






i 




p* _ 





where we have used the boundary condition 



p* 



(25) 



(26) 



(27) 
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This boundary condition expresses the basic assumption that 
triggering is attained when the sequence has been completed. To 
render the equations dimensionless, we introduce the following 
parameters: 



A 7 K* 



X-i 



X-2 



(28) 



The scaled (dimensionless) TCR triggering rate is then given by 
the following expression: 



where 



w = a ( 6 P<J + ? y off P*) 



w 



5 



1 + x + r + xr/ytin 
xr 

Ykin(l + x+ r) + xr 



(29) 

(30) 
(31) 
(32) 



The scaled TCR triggering rate w depends on ten dimensionless 
parameters (Table 2). 

STATISTICAL FORMULATION OF TCR DEGENERACY 

To express TCR degeneracy mathematically, we consider the dis- 
tribution of the triggering rate over the set of peptide ligands. 
This is just the set of Wy-values for a given TCR clonotype i over 
pMHCI ligands The distribution can be represented by plotting 
the probability P(w« > oo) as a function of w. Such a graph shows 
how many randomly selected peptides would have triggering rate 
w larger than a set value w. Let Tjj denote the mean dwell time 
of the TCR/pMHCI interaction for TCR clonotype i and pMHCI 
species j so that T^ 1 is the TCR/pMHCI off-rate (we have 
thus far suppressed subscripts for clone and ligand j to keep 
notation uncluttered). Arrhenius theory furnishes the following 
relationship with the dissociation energy barrier A Uif. 



Tij = T 0 exp(AU,j), 



(33) 



where To is the frequency factor and A Uy is expressed in Boltz- 
mann units. We assume that the dissociation energy barrier arises 
as a result of a large number of individual reaction steps at the 
TCR/pMHCI interface. If these combine additively, then the Cen- 
tral Limit Theorem implies that A Uy has a Gaussian distribution. 
Letting Uy = A Uy — ln{TR/To} we have ti y ~ N(— |x, a 2 ), where 
|x > 0 and a are the underlying parameters of the normal distribu- 
tion N. The assumption that the mean is negative is a consequence 
of the fact that for the vast majority of TCR/pMHCI complexes 
the mean dwell time is less than the typical time required to com- 
plete the chain of ITAM phosphorylations. It now follows that 
the dimensionless parameter a, associated with the TCR/pMHCI 
off-rate, is log-normally distributed. 

RESULTS 

We investigated the effect of variations of the total CD8 density 
on the functional sensitivity of hypothetical ligands with various 



Table 2 | Dimensionless (scaled) parameters that govern functional 
sensitivity. 

mj Scaled total pMHCI density 

Xj Scaled total CD8 density 

rj Scaled total TCR density 

a ScaledTCR/pMHCI off-rate without CD8 bound 

& Scaled pMHCI/CD8 off-rate with TCR bound 

v Scaled kinetic effect of pMHCI/CD8 interactions with and without TCR 

k Ratio of dissociation constants /Ci and K3 

y 0 ff Factor by which CD8 modulatesTCR/pMHCI off-rate 

Ykin Factor by which CD8 modulates theTCR/pMHCI affinity 

yfl Factor by which CD8 modulates the TCR triggering threshold 



TCR/pMHCI off-rates. All variables and parameters are dimen- 
sionless (scaled) in the model simulations. The scaled parameters 
are summarized in Table 2 and their scaling is defined in equations 
(14), (15), (17), and (28). 

LIGAND FOCUSING AND CD8-MEDIATED CONTROL OF DEGENERACY 

Figure 2A shows the scaled functional sensitivity w as a function 
of scaled total CD8 density xj with different scaled TCR/pMHCI 
off- rates a. A striking feature is the dip at xj ~ 50, suggesting that 
at this CD8 density the T-cell is minimally responsive. As the scaled 
total CD8 density increases above xj ~ 50, the scaled functional 
sensitivity increases for all hypothetical ligands. 

It can be observed that a ligand with a = 2.5 (solid line) 
and a ligand with a = 0.5 (dotted line) show opposing changes 
in the scaled functional sensitivity: a ligand that is less potent 
at low CD8 becomes more potent at high CD8 and vice versa. 
Hence, changes in CD8 expression levels can differentially affect 
the potency of ligands, each of which is potentially a strong ago- 
nist. In effect, the T-cell can tune in on a specific ligand and thus 
control ligand promiscuity. We call this the principle of ligand 
focusing. 

The corresponding degeneracy curves P(w > w) for different 
scaled total CD8 densities are shown in Figure 2B. The effect of 
the increase in co-receptor density on the sensitivity and degener- 
acy is indicated: sensitivity is a change in the horizontal direction 
(modulating the triggering threshold) and degeneracy is a change 
in the vertical direction (controlling cross-reactivity) . As the scaled 
total CD8 density increases up to %r~ 50, the degeneracy curves 
P(w > a>) move to the left from the curve without CD8 {xj = 0): 
degeneracy and sensitivity decrease. By contrast, as the scaled total 
CD8 density increases above xj ~ 50, the degeneracy curves move 
to the right, which means that degeneracy and sensitivity both 
increase. The overall effect is that the T-cell becomes more degen- 
erate as CD8 levels increase. Thus, in addition to the focusing 
effect, the co-receptor also governs the overall degeneracy of the 
T-cell. 

A high degree of degeneracy can increase the risk autoimmune 
disease. On the other hand, too low a degree of degeneracy could 
compromise the immune system's ability to mount a timely and 
efficient response. To analyze these risks, suppose that the T-cell 
is activated if its integrated TCR triggering rate exceeds a certain 
value, termed cellular activation threshold. A simple model is to 
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FIGURE 2 | (A) Scaled functional sensitivity w as a function of scaled total 
CD8 density x T for various scaled TCR/pMHCI off-rates a. (B) Degeneracy 
curves P(w > co) for various scaled total CD8 density x T . (C)The probability 
¥(w > to) as a function of CD8 density x T , at a set value of functional 



assume that the T-cell is activated if 

ZjT lWlj > W aa (34) 

where Zt is the presentation level of ligand j, 7/ is the duration of 
the T-cell:APC interaction, and W act is the activation threshold. 
For a given Zj and Tj, there is a critical w,j which is the mini- 
mum value required to satisfy equation (34). Suppose for instance 
that this corresponds to Wij = 0.12. The corresponding probabil- 
ities P(w> 0.12) are plotted in Figure 2C. Given the estimates 
of the TCR repertoire size, normal immune function is proba- 
bly confined to an operating range of probabilities 10~ 8 to 10 -5 . 
Figure 2C shows how the level of CD8 can regulate the respon- 
siveness to remain within this band; if the activation probability 
drops below this range, the risk of not responding to a pathogen 
looms, whereas at much elevated activation probabilities, the risk 
of autoimmunity is heightened. 

MODES OF CD8-MEDIATED MODULATION OF FUNCTIONAL SENSITIVITY 

Figures 3 and 4 exhibit the co-receptor effect on functional sen- 
sitivity in two parameter scenarios, where we assume a weak 



sensitivity co = 0.12. The operating range of the probability P is shown as a 
function of x T with dashed lines at P(w> 0.12) = 10"° and P(w > 0.12) = 10- 5 . 
Parameter values: 8 = 300, v = 0.5, n = 100, y lm = 0.5, y off = 0.5, y„ = 0.3, 
k = 5, m T = 10, r T = 10. The log-normal distribution has mean 2 and SD 0.2. 



(v = 0.05) or a strong (v = 300) kinetic effect of pMHCI/CD8 
interactions, respectively. When the kinetic parameter v is small, 
increasing the levels of CD8 on the T-cell surface leads to enhanced 
functional sensitivity of ligands with a > 1, as shown in Figure 3 A. 
By contrast, ligands with low off-rates (a < 1) become less potent 
when CD8 levels are increased. These opposite effects demonstrate 
CD8-mediated focusing on particular ligands. 

Figure 3B shows the effect of increasing CD8 levels on a lig- 
and that is optimal in the absence of CD8. When the modulatory 
effect of CD8 on the TCR triggering threshold is modest (y« 
is near 1), the ligand becomes less potent as levels of CD8 are 
increased. On the other hand, when CD8 has a strong effect on 
the TCR triggering threshold (low y#, since the time required to 
trigger the TCR/CD3 complex is shortened), increasing the levels 
of CD8 on the cell surface enhances the functional sensitivity to 
the ligand. 

The effect of CD8 levels on functional sensitivity is shown in 
Figure 4A. For each ligand, there is an optimal CD8 level which 
depends on that ligand's TCR/pMHCI off-rate a. This shows that 
the T-cell can favor the signaling strength of a given ligand by 



Frontiers in Immunology |T Cell Biology 



October 2013 | Volume 4 | Article 329 | 6 



Szomolay et al. 



Modulation of TCR degeneracy and sensitivity by CD8 




FIGURE 3 | Scaled functional sensitivity w as a function of scaled 
total CD8 density x T . (A) Curves for various values of the scaled 
TCR/pMHCI off-rate a. (B) Curves for various values of the factor y H by 



which CD8 modulates theTCR triggering threshold. Parameter values in 
(A,B): 8 = 2.5, v = 0.05, n = 100, y lln = 0.5, y off = 0.5, y„ = 0.5, k = 1 , 
m T = 10, r T = 10, oi = 1. 
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FIGURE 4 | Scaled functional sensitivity w as a function of scaled total 
CD8 density x T (A) Curves for various values of the scaledTCR/pMHCI 
off-rate a. (B) Curves for various values of the dimensionless factor y off 



expressing the modulatory effect of CD8 on theTCR/pMHCI off-rate. 
Parameter values in (A,B) 8 = 3, v = 300, n = 100, y km = 0.05, y olf = 0.2, 
y „ = 0.5, k = 2, m T =10, /y=10, a = 2. 



adjusting the density of CD8 molecules on its surface. As one 
would expect, the modulatory effect is most profound when the 
relative effect that the co-receptor exerts on the TCR/pMHCI off- 
rate is greatest; this is shown in Figure 4B. When CD8 alters the 
off-rate strongly, the functional sensitivity is strongly depressed 
when CD8 levels are increased beyond the optimal level, whereas 
for a moderate value of the modulatory multiplier y a g, the func- 
tional sensitivity remains at near-optimal levels when CD8 levels 
are increased. The value of y 0 S ma J b e expected to be different 
for different TCR/ligand combinations. In particular, when CD8 
makes a substantial contribution to the binding energy, the multi- 
plier Y 0 ff will be low, and the co-receptor role in governing ligand 
optimality will be more pronounced. 



CO-RECEPTOR MODULATION OF DEGENERACY VIA ALTERED 
pMHCI/CD8 BINDING AFFINITY 

Wooldridge et al. (5, 30, 31) have shown that (i) increased 
pMHCI/CD8 interaction results in enhanced recognition of 
pMHCI by cytotoxic T-cells and (ii) increased pMHCI/CD8 inter- 
action impairs pMHCI recognition specificity, suggesting that the 
pMHCI-CD8 interaction is essential in regulating the balance 
between optimal T-cell cross-reactivity and T-cell antigen speci- 
ficity (32). These findings suggest that an optimal pMHCI/CD8 
strength exists that yields maximum pMHCI sensitivity without 
loss of specificity. Motivated by these results, we consider a hypo- 
thetical scenario in which pMHCI mutant molecules with altered 
binding affinity for CD8 modulate TCR degeneracy. 



www.frontiersin.org 



October 201 3 | Volume 4 | Article 329 | 7 



Szomolay et al. 



Modulation of TCR degeneracy and sensitivity by CD8 



In keeping with (17), the following pMHCI mutants are 
considered: A245V representing weak pMHCI/CD8 affinity 
(JC D = 498u.M), wild-type (Ko = 137.1 u,M), Q115E represent- 
ing slightly enhanced affinity (Kd = 97.94 u,M), and A2/a3k b 
with enhanced affinity (Ko = 10.87 |xM). These values are based 
on Surface Plasmon Resonance experiments and should be 
regarded as "three-dimensional," relating to the ligands in solu- 
tion. However, the TCR/pMHCI interaction takes place in the 
"two-dimensional" environment of the T-cell:contact area, which 
essentially reduces spatial degrees of freedom of molecular motion 
and, moreover, introduces dynamics related to the forces that 
constrain the molecules to this environment, such as rotations 
with respect to the membrane plane, membrane fluctuations, and 
the translational motion of the membranes themselves (33). Fur- 
thermore, cooperativity interactions, such as the involvement of 
the co-receptor, may be profoundly altered. As a result, the rate 
constants can be markedly different in the "two-dimensional" 
environment; in particular, two-dimensional dissociation rates can 
be substantially faster (28). The ratio between the three- and two- 
dimensional affinities is a length measure, denoted h and called 
the confinement length (34, 35). Wu et al. (33) demonstrated that 
h is proportional to the range of motion available to the free forms 
of the interacting ligands along the spatial axis perpendicular to 
the two parallel membranes. Thus the inter-membrane separation 
distance provides an upper bound, and if the ranges of motion are 
broadly comparable, we may assume that the confinement length 
is roughly the same for all mutants involved. 

The two-dimensional dissociation constant for pMHCI/CD8 
interaction without TCR bound, K}, appears in the scaled para- 
meters k, xt, and v. Hence, by taking K3 to be proportional to 
Kd for a given pMHCI mutant, we can simulate the impact of 
the altered pMHC/CD8 binding affinity on TCR degeneracy and 
sensitivity. Figure 5A shows the scaled functional sensitivity w as a 
function of scaled TCR/pMHCI off-rate a for the four hypothetical 
pMHCI mutants. 

The co-receptor CD8 can modulate the specificity of antigen 
recognition, as shown in Figure 5B. Each TCR degeneracy curve 
corresponds to a given pMHCI mutant, where the one with the 
strongest pMHCI/CD8 binding affinity (A2/a3k b ) is the most 
degenerate, with the largest antigen sensitivity. The three regions 
are a schematic representation of the overall pattern of CD8 + T- 
cell antigen specificity, as defined by Cole et al. (32). With increas- 
ing strength of pMHCI/CD8 affinity, as indicated by the arrow, 
the recognition efficiency of partially CD8-dependent ligands is 
enhanced and the spectrum of CD8 + T-cell antigen degeneracy 
becomes wider. Enhancing the kinetic effect of pMHCI/CD8 inter- 
actions (setting v 3> 1) results in the reversed pattern, as shown in 
Figure 5C. 

Whereas in the MHC-limited kinetic regime, the behavior is as 
shown in Figure 5B, the degeneracy curves for the various mutants 
overlap in the TCR-limited regime. This suggests that excess levels 
of ligand, relative to the available levels of TCR molecules, can 
diminish the importance of the interaction between pMHCI and 
CD8. In principle, this endows the APC with a means to "override" 
the ligand focusing exerted by the T-cell, allowing a professional 
APC, such as a dendritic cell, to force a naive T-cell, with which it 
has conjugated, to be maximally degenerate. 



DISCUSSION 

The co-receptors CD4 and CD8 are glycoproteins that modulate 
the interactions of the TCR with pMHCI and pMHCII mole- 
cules, by binding to invariant sites on these molecules (10). It 
is well-established that the co-receptors differentially regulate the 
responsiveness of the TCR to the ligand and thereby modulate TCR 
specificity (1). In particular, CD8 is known to affect both the on- 
rate and the off-rate of the TCR/pMHCI interaction (17, 36). This 
allows the co-receptor to differentially regulate the strengths of 
the various potentially strong agonists of the TCR. This accords 
with the finding that the strength of pMHCI/CD8 interaction 
is a determinant of T-cell degeneracy (5). This ligand focusing 
effect remains to be observed experimentally, to the best of our 
knowledge. Perhaps this is only to be expected inasmuch as the 
experimenter has to search for ligands that are sub-optimal under 
standard conditions but become better or worse agonists when 
CD8 levels are manipulated. Research is presently underway to 
identify such ligands and we anticipate that the phenomenon 
will be confirmed and eventually emerge as a pervasive "design 
principle" of cellular adaptive immunity. 

Co-receptor-directed ligand focusing may allow the T-cell 
response to an antigen challenge to undergo an adaptive evolu- 
tion that would be functionally analogous to affinity maturation 
in B-cell immunity. Moreover, CD8 modulation could allow for 
an elevated degeneracy among the earliest responding clones. This 
would ensure that at least one or more responding clones are acti- 
vated sufficiently early in the course of the infection. Moreover, a 
gradual restriction of the degeneracy, coupled with an increase in 
functional sensitivity to the salient epitope, would then reduce the 
degeneracy of the response, which would gradually evolve from 
oligoclonal to one that is dominated by an optimally tuned single 
clone. 

Disrupting the pMHCI/CD8 interaction impairs the ability of 
T-cells to recognize antigens. In particular, T-cell activation can be 
abrogated if the pMHCI/CD8 interaction is blocked (32), whereas 
increases in pMHCI/CD8 affinity have the opposite effect (37, 38). 
The contribution of CD8 to increase functional sensitivity appears 
to be crucial for weaker agonists (36). A comprehensive evaluation 
of clonal CD8+ T-cell degeneracy using combinatorial peptide 
libraries and APCs expressing mutant HLA A* 0201 molecules 
with altered pMHCI/CD8 affinity has shown that the co-receptor 
enhances T-cell degeneracy by increasing the range of agonist lig- 
ands that can elicit T-cell activation (5). Furthermore, increasing 
the affinity of CD 8 for HLA A* 0201 by at least an order of mag- 
nitude resulted in the loss of cognate antigen specificity (5, 31). 
The affinity of the pMHCI/CD8 interaction may be directly linked 
to TCR degeneracy: increased pMHCI/CD8 interaction enhances 
CD8+ T-cell antigen sensitivity, but reduces CD8+ T-cell antigen 
specificity (32). This agrees with our main finding that variation 
of the co-receptor effect regulates the degree of T-cell degeneracy 
and antigen specificity. 

A cornerstone of the present theory is that a certain amount 
of degeneracy is unavoidable, in view of the vast universe of 
possible peptides and the relatively modest number of TCR clono- 
types that even a large mammal might be able to maintain in its 
standing repertoire. Moreover, salient epitopes, those associated 
with a disease state, and non-salient ones, such as self-peptides 
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FIGURE 5 | (A) Scaled functional sensitivity w as a function of scaled 
TCR/pMHCI off-rate ct. (B) Degeneracy curves P(w>m) for HLA A*0201 
mutants with altered binding affinity for CD8: A245V (dotted line), wild-type 
(semi-dashed line), Q115E (dashed line), andA2/a3k b (solid line). The three 
regions represent the overall pattern of CD8+ T-cel! antigen specificity and the 



for which immune tolerance is required, will of necessity be 
"finely interleaved" subsets of the peptide universe (a mathe- 
matician would say that one subset is "dense" in the other), 
lest the tolerant subset forms a target for the rapidly evolving 
pathogens: the system cannot work if molecular mimicry is read- 
ily attained. From these two premises, it follows that a TCR must 
be degenerate, and also that this degeneracy must be suscepti- 
ble to exquisite modulation. Against this line of reasoning a case 
could be made that the size of the ligand universe is effectively 
much smaller. For instance, if one considers n-mer peptides that 
are anchored to the MHC binding groove at a positions, and 
the region of the TCR that interacts with the peptide (roughly 
speaking, the CDR3 loop) makes contacts with c of the amino 
acid residues in the m-mer (so that 0 < m < n — a), then there 
are 20 m (n— a)\l{m\{n— a— m)\) effectively distinct pMHCI lig- 
ands as seen by the TCR. To give an extreme example, with n = 9, 



arrow indicates the strength of pMHCI/CD8 interaction. (C) Degeneracy 
curves V(w > to) for HLA A*0201 mutants with altered binding affinity for CD8. 
The parameter values are the same as in A except for v= 10. Parameter 
values in (A,B): 5 = 0.2, v = 0.05, n= 100, y lin = 0.5, y„„ = 0.5, y B = 0.2, k = 1, 
m T = 10, r T = 10, x T = 10. The log-normal distribution has mean 5 and SD 0.5. 



a = 2, and m=l, this works out as just 20 x 7 = 140 distinct lig- 
ands. Perhaps the estimate m = 4 is more realistic: this gives only 
20 4 x 35 = 5.6 x 10 6 functionally distinct ligands, which is of the 
same order as the TCR repertoire size. Whereas there may be 
some merit to this argument, its underlying image, essentially of 
CDR3 as a tape recorder head that interacts with only m amino 
acid residues and is indifferent to the n — m other ones, is a 
gross oversimplification. The physical behavior of the m amino 
acid residues at the contact sites cannot fail to be influenced 
by the n—m remaining ones (including the a anchor residues). 
This is overwhelmingly apparent not just from the basic prin- 
ciples of molecular dynamics, but also from the typical results 
obtained with combinatorial peptide library scans. A case in 
point is the finding that changes at the anchor position result in 
changes in the center of the peptide and therefore influence TCR 
binding (39, 40). 
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The present model indicates that intermediate levels of CD8 are 
associated with the lowest functional sensitivity. This suggests the 
following mechanism to maintain quiescent (naive) T-cells in a rel- 
atively unresponsive state. When the T-cell receives the appropriate 
stimuli, it either up-regulates or down-regulates the co-receptor 
and a specific subset of its potential agonists "comes into focus." 
Such signals are known to be transmitted via cytokine profiles in 
the T-cell's surroundings (15) as well as costimulatory receptor- 
mediated signals transmitted by professional APCs (41). When no 
harm is detected, the default response of the naive T-cells would 
be to "de-tune" whenever a strong signal is registered. Detuning 
of T-cells via alterations of CD8 expression levels, under control 
of cytokine stimuli, has been reported (42) and the connection 
between functional sensitivity, tolerance, and CD8 expression lev- 
els is well-established ( 1 4, 43, 44) . On the other hand, when harm is 
detected and transmitted via a pro-inflammatory cytokine profile, 
the T-cell's response would invert and the tendency would become 
to "tune in" to any supra-threshold stimulation. On a molecular- 
cellular level, this would involve a scanning mechanism whereby 
the stimulus would make the cell enter a mode in which it gradu- 
ally alters the CD8 expression level whilst the received TCR signal 
(which will shift in magnitude as the CD8 level changes) feeds 
back onto this pathway. To the best of our knowledge, the molecu- 
lar details of such a regulatory pathway have not been elucidated to 
date. However, we believe that it is well within the regulatory capa- 
bilities of cellular signaling networks; we have previously discussed 
similar mechanisms in more depth (45, 46). 

Whilst the model includes the key components of TCR trig- 
gering, many important aspects have been omitted. In particular, 
we have neglected the spatial dynamics of TCR, CD8, and pMHCI 
within the immunological synapse, where the relative concentra- 
tions of p56 lck and CD45 will determine how quickly partially 



phosphorylated TCR/CD3 will reset to the basic state (47). A 
functional consequence of the exclusion of phosphorylases, for 
instance, could be the ability of the TCR/CD3 complex to be trig- 
gered over subsequent interactions with ligands that would in a 
normal context only have weak functional sensitivity. 

In addition to kinetics of the interaction between TCR, MHCI, 
and CD8, we have only considered ITAM phosphorylation steps. 
It is well known that phosphorylated ITAMs orchestrate the acti- 
vation of the Src-related protein tyrosine kinases which initiate 
TCR signaling. These kinases induce tyrosine phosphorylation of 
several polypeptides, including the transmembrane adaptors. Pro- 
tein tyrosine phosphorylation subsequently leads to the activation 
of multiple pathways such as ERK, NF-kB, and NFAT (48, 49). 
Moreover, negative regulation of TCR signaling is key to avoiding 
hyper-activation. Notwithstanding the additional layers of com- 
plexity which our simple model ignores, we believe that the system 
of two linked proofreading chains as presented here does capture, 
in a qualitative sense, the essence of TCR triggering. 

In summary, the present findings suggest that the co-receptor 
CD8 can differentially modulate functional sensitivity to its poten- 
tial agonists, thereby modulating TCR degeneracy in a tunable 
fashion. The ligand focusing mechanism would allow each T-cell 
to have a wide range of potential agonists, even while only one 
of these would be a ligand of high functional sensitivity at any 
particular moment in time. 
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